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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

CORRELATION OF PREDICTED AND EXPERIMENTAL LATERAL 
OSCILLATION . CHARACTERISTICS FOR 
SEVERAL AIRPLANES 

By Donovan R. Heinle and Waiter E. McNeill 

SUMMARY 


An investigation of the accuracy of a presently used method for 
predicting the dynamic lateral oscillatory characteristics of airplanes 
has been made. This study involved the comparison of contr ole-fixed 
flight-test data from several current fighter-type jet airplanes with 
the calculated values of period and damping for these airplanes. These 
calculations were made from geometric and mass characteristics and aero- 
dynamic derivatives derived from geometric characteristics and wind- 
tunnel data. 

For the airplanes included in this investigation, it was found that 
the dynamic lateral oscillatory characteristics could be predicted with 
reasonable accuracy by 1186 of the geometric and mass characteristics and 
static wind-tunnel data. Where discrepancies did exist between experi- 
mental and calculated data, they were attributable to control motion or 
to parameter inaccuracies in the calculations. It was found that, for 
airplanes involving deviations from the conventional configurations, 
special techniques for obtaining the dynamic derivatives such as dynamic 
wind-tunnel tests are advisable. 


INTRODUCTION 


f 


With the increasing complexity and cost of the modern airplane, it 
becomes desirable to predict as far as possible in advance of its com- 
p letion approximately what the flight stabili ty characteristic s win be . j 
"Tims any expected deficiencies may be noted &ncL xheir possible corrections 
may be applied early in the development and construction of the airplane \ 
to minimize time loss and expense. The extension of flight to the higher 
altitudes and speeds experienced in recent years has made it an increas- 
ingly greater problem to produce aircraft that perform within the limits 
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of desira ble-stability characteristics . so the prediction of thes 
stability characteri sticg _hecameB-^iore 1 Tnpor.t>aB&. 


:sej 


The purpose of this investigation was to find if a currently used 
method of calculation of the lateral oscillation characteristics of an 
airplane is sufficiently accurate. This study concentrated on the short 
period oscillation, commonly referred to as the "Dutch roll" mode, and 
did not include the investigation of snaking which might appear on some 
aircraft. The calculated values of period and time to damp to half 
amplitude of the lateral oscillations (controls held fixed) were compared 
with experimental flight-test values to see how closely these motions 
could be predicted. Any discrepancies which appeared between calculated 
and experimental data were examined for possible causes such as aero- 
elasticity, unknown Mach number or altitude effects, aerodynamic lag, or 
incomplete derivative estimation. 


It was decided to concentrate, in this Investigation, on the smaller 
fighter-type airplanes presently flying. At the beginning of the inves- 
tigation published data were not complete on most of the aircraft. 
Therefore, it became necessary to gather flight data from the manufac- 
turers of the individual airplanes. In addition, wind-tunnel, geometric, 
and mass data for the comparable configuration were requested. The 
flight data used were obtained from tests conducted by the Armed Services, 
by the manufacturer, and by the NACA. The wind-tunnel data with the 
inertia and geometric characteristics were analyzed by a more or less 
standard calculation procedure. 


The methods of calculation and the results of this survey are 
reported herein. 


DESCRIPTION OF AIRPLANES 


Geometric and mass characteristics and aerodynamic derivatives for 
each airplane studied are listed in table 1 for two altitudes and one 
Mach number. Two-view silhouettes are presented in figure 1. 

The airplanes studied have jet engines. Airplanes A, D, E, F, and 
H are twin engined, while the rest have single engines. As can be seen 
in figure 1, the nin e airplanes exhibit a variety of wing plan forms and 
fuselage shapes. Airplanes A, C, D, and E have straight wing and tail 
surfaces. The vertical stabilizer of airplane C has a swept-back leading 
edge. Airplane B has straight surfaces except for the horizontal tail, 
which is swept back 35° • Airplanes F, G, and IT halve all surfaces swept 
back approximately 35° • Airplane H is a semitailless design having no 
horizontal tail surfaces and with vertical surfaces on the wing. Air- 
plane I has an inverse taper, wing swept back 40° at the midchord line. 
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Its horizontal and vertical tail surfaces are swept back 40° and 33° > 
respectively. Airplanes A, B, and G have low wings; airplanes C and F 
have low midwings; airplanes D, H, and I have midwingB; and airplane E 
has a high midwing. Six of the airplanes (A, C, F, G, H, and I) are 
single-seat fighters; B, D, and E are two-seat night or all-weather 
fighters. 


PROCEDURE 
Data Compilation 


The flight data available from reports by the Armed Services and 
NACA were not sufficient to conduct an adequate survey. Consequently 
the manufacturers of each airplane were contacted to obtain the neces- 
sary information. The manufacturers were requested to furnish flight 
data of the type desired and information on the mass and geometric char- 
acteristics and wind-tunnel tests comparable to those of the test vehicle. 
The wind-tunnel test results were available in reports (some unpublished) 
from the NACA and other wind-tunnel facilities. The desired flight-test 
information was requested to be in the form of time histories of controls- 
fixed lateral oscillations at several altitudes over a range of speeds. 
These time histories were to furnish, more specifically, the period and 
time to damp to half amplitude of the lateral oscillations. 

The manufacturers furnished flight data together with static wind- 
tunnel data and geometric characteristics for airplanes A, C, D, E, F, 
and H. Additional flight data were obtained from Armed Services reports 
for airplanes C and D. Wind-tunnel data and geometric characteristics 
were furnished for airplanes B, G, and I by the manufacturers. Flight 
data for airplane B came from Armed Services reports and for airplanes 
G and I from NACA tests. 


Calculations 


t 

Appendix A presents a complete list of the symbols used in the cal- 
culations and the presentation of the data. Figure 2 shows the system 
of axes used for the calculations. The method of calculating the period, 
P, and the time to damp to half amplitude, Tj y 2 , is outlined in detail in 
appendix B. References 1 through k- were used to estimate the rotary 
stability derivatives, neglecting frequency effects (reference 5 ) . 
Reference 6 was used to estimate the effect of Mach number on the rotary 
derivatives and the low-speed static stability derivatives. The airplane 
was assumed to be rigid (no aeroelastic effects), and all controls were 
assumed fixed. 
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For additional methods of obtaining the lateral stability deriva- 
tives, see reference 7* 


RESULTS 


Comparison of the experimental and predicted values of period and 
time to danrp to half amplitude is shown- for each airplane in figures 3 
through 11. The flight- test data are shown as individual points and 
calculated data are presented as curves. The altitudes range from 
8,000 feet to 35*000 feet. It will be noted that some flight-test point 
symbols are filled indicating controls-fixed data. This was done only 
when it was definitely known by records or other means that the control 
surfaces were fixed. Since - controls-fixed data had been requested for 
this investigation, they were assumed controls-fixed (indicated by the 
half— filled symbols) unless definitely known to be otherwise, as proven 
by control position records or by the method of test. If control move- 
ment was indicated, the symbols are open . 

Most of the flight-test data available from the manufacturers had 
been gathered incidentally to other tests rather than from testing spe- 
cifically for controls-fixed lateral-oscillation characteristics. As a 
result, the data available for some of the airplanes were rather incom- 
plete . 

No effort was made in this investigation to study snaking, a char- 
acteristic indicated by flight records of some of the airplanes at high 
altitudes. If the initial amplitude of the Dutch roll oscillation was 
large enough, the value of Ti/z could be obtained and would be essen- 
tially Independent of the snaking. Otherwise, only the period was 
recorded. 


DISCUSSION 


Based on the present and previous lateral stability, investigations, 
it appears that ability to predict controls-fixed period^ P, and damping, 
l/T 1/s , to within about +20 percent is sufficiently accurate to be use- ' 
ful in the design of tactical aircraft. Consequently, this figure will 
be used to distinguish between satisfactory and unsatisfactory results 
in the following presentation. . ... 

Several factors determined this choice. It is doubtful that current 
or future flying-qualities requirements involving lateral-oscillation 
period and damping can define satisfactory and unsatisfactory behavior 
to any greater accuracy than ±20 percent. Experience has indicated 
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(see reference 8, for example) that pilots generally do not discern or 
react to changes in period and damp in g of much less than ±20 percent. 
Scatter in flight data can be attributed to small variations in speed 
and altitude, recording inaccuracies, and errors or generalizations 
involved in data reduction. It has been found that this scatter can be 
held below ±20 percent only in a carefully performed experiment. 

In figures 12 and 13 are plotted the values of measured period and 
damping (reciprocal of time to damp to half amplitude) against the pre- 
dicted values for each of the airplanes. Also shown on these plots are 
the lines of perfect agreement and ±20 percent disagreement. Generally, 
it may be seen in figure 12 that the period could be predicted with 
reasonable accuracy. Xn figure 13, it may be noted that with the excep- 
tion of airplanes B, E, and I, the prediction of damping was satisfactory. 
The reasons for the poor correlation for these three airplanes will be 
discussed later. 


Airplane A 


It may be seen in figure 3 that for this airplane the correlation 
between experimental and predicted period at 10,000 feet altitude and 
time to damp to half amplitude at 10,000 and 30*000 feet is good. The 
comparison of experimental and predicted period at 30,000 feet altitude, 
however, shows poor agreement. A general increase of the predicted curve 
by 30 percent would bring about good correlation. It is not known what 
caused this disparity between the measured and predicted data at 30,000 
feet altitude. Comparison with the data of the other airplanes included 
herein of similar configuration shows that the measured period of air- 
plane A increases by 70 to 80 percent in going from 10,000 to 30,000 feet 
altitude; while the measured values of period for the other airplanes 
increase by only 40 to 50 percent. This fact might cause doubt of the 
accuracy of the measured data except that the measured data available 
appear to be repeatable. There is the possibility that a moving free 
rudder could be the cause of the increase in period. 


Airplane B 


The data available for airplane B were too sparse to provide a ready 
means of comparison between the predicted and measured values of period 
and damping. Figures 12 and 13 indicate good agreement among the values 
of period but poor correlation among the damping values. It will also be 
noted (see fig. 4) that the damping decreases markedly in the low Mach 
number range . 
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At the time the calculations were made, wind-tunnel sideslip data 
were not available for the complete airplane configuration tested in 
flight. An increase of wing geometric dihedral to 7*0° for airplane B, 
as compared with 3.6° for the preceding airplane of the series (for which 
low-speed wind-tunnel sideslip data wera-.available) , accompanied by a 
reduction in wing thickness, required that corrections be made to dihedral 
effect, Cjg, The low-speed values' of Cnp and were taken directly 

from the wind-tunnel data available and all rotary derivatives were com- 
puted by the methods of appendix B . 

The high dihedral, effect of this airplane is the main cause of the 
low damping present and is partly responsible for the poor correlation 
between the predicted and measured damping. 


Airplane C 


The correlation for the period and the damping is good at the 
10,000-foot altitude (fig. 5 ). Predicted and experimental values of 
damping do not agree’ well at 30,000 feet. The measured values show 
better damping than predicted. The scatter in the data is a little 
larger than usual. Which could cause some of the discrepancy. 


Airplane D 


The data available for airplane D indicated good correlation (see 
figs. 6, 12, and 13)* The test data for damping at 30,000 feet (not 
shown) indicated considerable scatter and were considered unreliable for 
comparison. Most of the flight-test data were obtained with the stick 
held fixed but with the rudder free. Small oscillations of the rudder 
which were present could possibly be responsible for the slightly higher 
experimental times to damp to half amplitude as compared with predicted 
values. The necessity of actually having the rudder . looked during flight 
tests cannot be overstressed since even small oscillations in proper 
phase allowed by cable stretch from locked controls can change the damp- 
ing considerably (see reference 9> for example). 


Airplane E 


The correlation of predicted and measured period of airplane E is 
generally satisfactory (figs. J, 12, and 13)* The damping correlation 
is not satisfactory. Figure 7 shows variations in the experimental data 
which accompanied the noted changes in conf iguration during the 
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flight-test program. This airplane, as tested, was subject to snaking 
oscillations so that some of. the flight data for dampi ng were unsatis- 
factory for comparison. 

The desirability of obtaining more extensive dynamic wind-tunnel 
data was demonstrated in the development testing of airplane E. When 
the early flight tests indicated that the damping was lower than pre- 
dicted, a model was placed in the tunnel by the manufacturer for further 
tests. This model was spindle mounted, permitting it to be oscillated 
in yaw. Tests of this arrangement also gave less damping than predicted 
from static wind-tunnel tests. Wake surveys subsequently made in the 
region of the tail showed that the fuselage was acting as a lifting stir- 
face in yaw, with vortices from the fuselage trailing near the tail. 

The high position of the vertical tail, however, placed it only in the 
upper part of this vortex flow, which is destabilizing, that is, it forti- 
fies the side flow in the direction in which the tail is moving. The 
destabilizing effect of the vortices could be reduced by use of a lower 
vertical tail position, an expensive problem in redesign after an air- 
plane is built. 

The above change in damping would be more properly termed Cnp 
than a change in Cn r (see symbols, appendix A). In that case, the 
total damping-in-yaw derivative should include the effects of both 
C nr and Cng • In the calculations for airplane E, the results of which 
are shown in figure 7 , Cn r was modified to include the approximate 
effect of Cnp ns indicated in appendix B. This modification increased 
the predicted values of time to damp to half amplitude by 38 to 89 per- 
cent. Plotted in figure 7 are curves of time to damp to half amplitude 
with Cn r modified and Cn r unmodified (dashed line) for one altitude 
( 8,000 feet). 

It would appear then that care in estimating the rotary damping in 
yaw must be exercised if accurate values are to be computed, especially 
if the static wind-tunnel tests indicate the presence of large sidewash 
variations at the vertical tail. 


Airplane F 


The majority of the flight-test points available for airplane F 
(fig. 8 ) were at 10,000 feet altitude. These indicate a slightly shorter 
period and slightly longer time to damp to half amplitude than the cal- 
culated values. Generally, there is good agreement between predicted 
and measured values of period and damping for the three altitudes shown 
(see figs. 12 and 13 ) . Some of the flight-test points for damping were 
not plotted due to incomplete or too short time histories available. 
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Airplane G 


Airplane G allows very good correlation between predicted and meas- 
ured values of period (see fig. 12). Generally, the correlation of 
damping is good (fig. 13) with a few points falling outside the ±20 per- 
cent boundaries. It may be seen in figure 9 that most of these points 
are at the limits of the scatter present. A more complete description 
of the correlation between predicted and measured lateral stability 
characteristics for airplane G may be found in reference 10. 


Airplane H 


The six experimental test points available for airplane H (fig. 10) 
indicate fairly good agreement between the predicted and measured values 
of period and damping. These few flight-test points, however, are too 
sparse to provide a definite comparison. 


Airplane I 


Reference to figure 11 shows, for aipplane I, marginal agreement 
between measured and calculated values of period. The prediction of 
time to damp to half amplitude is poor. Predicted values of the period 
are somewhat lower than the flight results over the entire Mach number 
range for the altitudes of 10,000, 20,000, ' and 30,000 feet, with a maxi- 
mum disagreement of —20 percent.. The predicted variation of T x / 2 with 
Mach number is generally in poor agreement with experimental results, 
and shows an over-all decrease with increasing Mach number much greater 
than is shown from flight tests. At 20,000 and 30»000 feet, the agree- 
ment between measured and calculated damping is very poor, particularly 
in the low Mach number range tested at each altitude. At 30,000 feet, 
for example, at a Mach number of 0.55; the calculated value of T 1( / 2 is 
approximately 100 percent greater than the experimental value. Because 
of the unusual configuration of this airplane (swept, inverse -tapered 
wing, and short tail length), it is likely that present methods of com- 
puting the rotary derivatives (Cn p ; Cn r , etc.) from static wind-tunnel 
data and geometric characteristics are not applicable. Period and time 
to damp to half amplitude were recomputed for the range of Mach numbers 
considered at 10,000. feet, neglecting Cn p . The results of these cal- 
culations are presented in figure'll and show good agreement between 
predicted and experimental damping at the lower Mach numbers, indicating 
that Cn p may be much smaller in magnitude "throughout the lift coeffi- 
cient range than was computed. Direct measurement of the rotary deriv- 
atives in curved- and rolling-flow stability tunnels would probably 
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furnish data more reliable, for use In calculating the damping of this 
airplane . 


CONCLUDING REMARKS 


Where a sufficient (enough to show repeatability) number of good 
flight-test points were available, the correlation of predicted and 
measured lateral-oscillation characteristics was usually satisfactory 
for the smaller fighter-type jet airplanes investigated. In all cases 
studied, the order of magnitude of the damping could be predicted. In 
general, the prediction accuracy with regard to the damping was much 
greater for well-damped airplanes, indicating the error may be a given 
uncertainty rather than a given percentage. 

A configuration with which the designer is familiar, such as one 
having only slight changes from previous airplanes, may have its lateral 
characteristics calculated fairly accurately by this method. However, 
a radical change in a model or a new configuration makes more extensive 
wind-tunnel tests highly desirable and In many cases necessary. These 
tests should include at least determination of Cn r and Cng as well as 
the usual static data. Good experimental values of the other rotary 
derivatives are also desirable. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif. 
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APPENDIX A 
SYMBOIS 


The stability system of axes used, in this investigation is defined 
as an orthogonal system having its origin at the center of gravity and 
in vhich the X axis extends along the flight path, the Z axis is in the 
plane of symmetry and perpendicular to the X axis, and the T axis is 
perpendicular to the plane of symmetry. See figure 2. 

b ■wing span, feet 

i w wing incidence, degrees 

l distance of center of pressure of vertical tail aft of the 
Z axis, feet 

m mass, slugs 

p rolling angular velocity about the X axis, radians per second 
q dynamic pressure ^ipV 2 ^ , pounds per square foot 

r yawing angular velocity about the Z axis, radians per second 

t time, seconds 

z distance of center of pressure of vertical tail from X axis, 
positive when center of pressure is above X axis, feet 

L rolling moment about X axis, foot-pounds 

M Mach number 

N yawing moment about Z axis, foot-pounds 

P period, seconds 

S wing area, square feet 

T^ 2 time for amplitude of oscillation to decrease by one-half, seconds 
V velocity, feet per second 

W weight, pounds 

T lateral force along T axis, pounds 
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a 

P 

7 

e 

*1 


P 

cr 

9 

t 




Cn 


C T 





angle of attack of fuselage reference line, degrees 
angle of sideslip, radians 

angle between flight path and horizontal, degrees 

angle "between reference axis and principal longitudinal axis, 
positive when tbs reference axis is above the principal axis 
at the nose, degrees 

inclination of the principal longitudinal axis of the airplane 
with respect to the flight path, positive when the principal 
axis is above the flight path at the nose, degrees 

relative density factor f 

\pSb 

mass density of air, slugs per cubic foot 
sidewash angle, positive for positive side force, radians 
angle of bank, radians 
angle of yaw, radians 

lift coefficient [ 

\ qS 

rolling-moment coefficient 
yawing-moment coefficient 
lateral-force coefficient 

ac z 
sp 

<*n 

dP 

dC n 

a(dp/dt) 

dCy 
dP 
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d(pb/2Y) 

ten 

5(pb/2Y) 

5Cy 

d(pfc/2Y) 

SC 2 

a(rb/2Y) 

Sc n 

d(rb/2Y) 

aGy 

a(rb/2Y) 

radius of gyration in roll about principal longitudinal axis, feet 
radius of gyration in yaw about principal vertical axis, feet 


K x nond i mens i onal radius of gyration about longitudinal stability 

axis («/ K Zq 2 cos 2 t) + K Zq 2 sin 2 t) ) 


Kj, nondimens ional radius of gyration about vertical stability 
axis (^/K Zo 2 cos 2 r[ + £x 0 2 sin 2 ^ .) 


nondimens i onal product— of— inertia parameter 

t(Kz 0 2 “ cos ti B±n ti 3 
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APPENDIX B 
CALCULATION METHODS 


In making the calculations of this investigation, the airplane vas 
treated as a rigid "body and values of period and damping "were found "by 
use of the following three linearized nandimens i anal equations of lat- 
eral motion referred to the stability system of axes: 


Boll 


(^-i I C Zp d) q) + (airf - ^ C 2r D ) C *f3 P =0 

T ear 

(^xzD 2 - I 2 Cn p 33 j q, + ^l z ^ a - 1 1 Cb, ^ P = ° 


Sideforce 


(- | Cy p D - 2 Cl) cp+ -|c Tr ^D — Jcl tan 7 * + (%iD - g C Yp ^) (3 = 0 


where 


D = 


d_ 

dt 


By use of these equations the stability quartic vas obtained. This 
was solved by synthetic division to obtain values of the period snd time 
to damp to half amplitude as follows: 


P 


2rt 

b 


^•i/ 2 


-O.693 

a 


where a «nd b are the real and Imaginary parts of the complex roots, 
respectively. The data to ascertain the coefficients of these equations 
were supplied partially by the manufacturers concerned and partially by 
normal methods of computation from basic data. The .weight and mass data 
were obtained from the man ufacturers together with values of the static 
lateral stability derivatives (C , Crip , Cy^ ) which had come from 

wind-tunnel tests of the model most closely representing the airplane 
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upon which flight teste were made. The rotary or dynamic derivatives 
were obtained by NACA methods of calculation (outlined herein) using the 
static derivatives and the airplane geometric characteristics as a base. 
The derivatives could be approximated by considering the geometric 
theory alone, but the results of wind-tunnel tests provide more accurate 
data to use. Kecent wind-tunnel tests using rotary balances or curved 
or rolling flow have been successful in providing empirical data useful 
for specific configurations. 

For these calculations no attempt was made to include a fuselage 
effect. It is usua l ly considered negligible for conventional airplanes . 
Therefore, the contributions of the wing and the vertical and horizontal 
tails to the various derivatives were obtained and then algebraical l y 
added. Generally, the equations of reference 1 modified for Mach number 
effects by the methods of reference 6 were found to give reasonable 
results for the wing contribution. 


Damping in Roll, 


The values of both the wing and horizontal tail contributions 
to CUp were obtained by the methods of reference 2. The charts pre- 
sented in that reference provided a rapid means of estimation for wide 
ranges of aspect ratios, taper ratios, and positive and negative sweep 
angles . 

The vertical tail contribution to C 2 p as well as Cj^ and Cy^ 
was estimated by the equations presented in reference 4 as follows: 

(f cos a ~i Blna )K» 008 a “ 1 61 “ a ) ~ ] 


In this equation, ckr/d(pb/2V) is termed the "sidevash factor." 
Values of this factor found to apply over a moderate range of angles of 
attack and for wings with and without sweep and a further discussion may 
he found in reference 3- Generally, it was found that the introduction 
of this sidevash factor served to reduce the vertical tail contribution 
so that Ci^ far the total airplane could be approximated by the value 


of 



for the wing alone. Values of Cy^ 


equations were obtained from wind-tunnel data. 


in this and the following 
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Yawing Moment Due to Rolling, Cjjp 

The contriTmtion of the wing to was determined by the methods 

of reference 1. These v alu es were further corrected for Mach number 
effects by equation (3) of reference 6. The horizontal tail contribution 
was considered negligible. The vertical tail contribution was determined 
by the following equation from reference 

% t ■ -° T P t (? 81,1 a + i 008 a ) [ 2 (t 008 a - i 8111 “)" sf^rj] 


Side Force Due to Rolling, Cyp 


The side force due to roll of the wing is dependent upon the sweep. 
Straight wings are assumed to have no, or negligible, side force due to 
roll. Swept wings do contribute a value to this derivative. Yalues for 
the wipg were obtained from the methods of references 1 and 6; the hori- 
zontal tail was neglected. The following equation of reference 4 was 
used to estimate Cyp due to the vertical tail: 


CYp t = Cy, 


Pt 


£ cos 
b 


a — 1 sin a 

D 


■) 


for 

3(pb/2Y) J 


In the calculations for several of the airplanes included herein, 
values of Cyp and Cy^. were not calculated or used since they have a 
negligible effect an the period and damping. 


Rolling Moment Due to Yawing, Cj r 


The horizontal tail was neglected in obtai ning the yawing dynamic 
derivatives. The contribution of the wing to Cj r was obtained by means 
of references 1 and 6. The vertical tall contribution was obtained by 
the following formula (reference 11): 

C 7 = -2 if z_ l sina \ c 
z r t b V b b J x Pt 
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Damping in Yaw, C^. 


No Mach, number correction was applied to the values of Cjjj. for 
the wing as obtained from reference 1. The vertical tail contribution 
to C nr was obtained by the following formula (reference 11): 


As noted in the discussion, the damping of airplane E was found to 
be less than that calculated by use of static wind-tunnel data. Conse- 
quently, for the study of this airplane, the calculated value of Cn r 
for the total airplane, at a given Mach number and C^, was reduced in 


magnitude by an amount equal to the difference between C 


n Ft 


as calcu- 


lated and Cn_ as obtains d--fn the wind-tunnel tests of the spindle— 


±r t 

mounted model of airplane E. The difference between C- 


nr t 


as calcu- 


lated and as measured in the wind tunnel is approximately equal to C -np 

These modified values of Gq^. were used to obtain the predicted curves 
of figure 7. 


Side Force Due to Yawing, Cy r 

The wing contribution to Cy r was obtained by the method of refer- 
ence 1 with no Mach number correction. Vertical tail values of Cy r 
(reference 11) were obtained by: 

* -* i 


It should be noted that the success of these calculations in pro- 
ducing good values of the rotary derivatives is dependent upon accurate 
static wind-tunnel data. 
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TABLE I.- PARAMETERS USED IK CALCULATIONS FOR AIRPLANES INVESTIGATED 
AT TWO ALTITUDES AND ONE MACH NUMBER 

[M » 0 . 50 ) 



A 


1 






X 




G 




I 


Altitude, ft 

10,000 

30,000 


20,000 



10,000 


0,000 


10,000 

30,000 

10,000 

33,a» 



10,000 

30,000 

V 



i5,ers 

i3,arfi 



16,100 


31,000 

31,000 

14,995 

14,995 

18,500 

12,300 


17,200 

19,535 

19,535 

b 

11.68 

41.62 

37- fi 3 

37-63 

35.25 

33.25 

50.00 

30.00 

5e.oo 


39.64 

39-66 

37.12 

37-12 

38, 6T 

36.67 

31.83 

31.83 

fl 

291.09 

294.09 

232.60 


250.00 


WQ.00 


606.00 

606.00 


390.00 

867.90 

tfT.90 

196.00 

196.00 

320.00 

320,00 

w/b 

P-.O 

k7.6 

65.6 

65.6 

52.0 

3B.0 

*5.8 

15.2 

31.1 

31.1 

10.3 

18.3 

43-4 

43.4 

36.1 

31.7 

61.0 

$1,0 

M 

156.0 

k3k.5 

171.9 

171.9 

tol.i 

lob.l 

562.1 

362.1 

963.63 

963.63 

526.3 

586.3 

388.5 

388.5 

567.3 

534.7 

607,24 

$07.81 

1* 

81.47k 

39.936 

13.303 


s6.ni 

51-581 

I6.005 

31.613 

16.368 

31.398 

21.672 

12.808 

20.702 

19.393 

17.143 

31.3?B 

34,603 

$0,319 

°L 


0,433 

0.256 

0,506 

0.201 

0.173 

0.178 

0.112 

0.186 

0.163 

0,107 

0.143 

0.171 

0,499 

0.15Q 

0.328 

0.848 

0.55$ 

a 

1.5a 

1.00 

1.37 

'2.61* 

2.11 

5.06 

-2.15 

0.35 

2.31 

$.21 

2.55 

3.91 

1.97 

6.31 

8.54 

5.43 

8.95 

8,15 

a, 

0 

0 

0 

« 

0 

0 

1.09 

1.09 

1.5 

1.3 

1.00 

1,00 

root + 1 
tip - 1 

root + 1 

up - 1 

0 

0 

0 

0 


-a-37 

0.03 

0.1k 

l.U 

-1.31 

1.11 

-1.30 

-1.70 

-2,76 

0.91 

-0.90 

2,16 

1‘MI 

4.01 

1.6? 

4.36 

l.k? 

6. $5 

« 

3.9? 

3.9? 

1.23 


3.95 

. 3.95 

2.05 

2,05 

5.30 

3.30 

3.45 

3.43 

2.50 

2,30 

g.ar 

ifKHTJ 

1.30 

1.50 


0.01»1 

0,015kT 

0,01003 

O.Q10U 

0.01529 

0.01326 

0.0UI6 

0.01132 

O.O3260 

0.03233 

0.W-32J? 

0,01336 

0.01333 

0.01367 

0.02801 

0.08620 

O.QGJOfl 

0,08583 


0,0*063 

0.01097 

S3 

0.07095 

0.01970 

0.01979 

0.03630 


0.05819 

0.03631 

0.07522 

0.07312 

0,04332 

0.043)8 

0.06635 

0.06503 

0.06666 

0.06731 


-0.00103 

0.0000! 

0.00013 

0.00130 


0.0006T 

-0.00213 

-O.OOOBl 

-0,00121 

O.OOOll 

-0.000*7 

0,00266 

-0.00026 

0.00806 

0.00112 


0.00161 

0.00733 

% 

-0.4498 

-0.4400 

-0.130 

-o.k 29 

-0.166 

-0.161 

-0.170 

-0.176 

-0.371 

-0.378 

-0,3539 

-0.3433 

-0.338 

-0.357 

-0,261 

-0.897 

-0.314 

-0.268 



0.1250 

0.307 

0.12J 

0.093 

0.351 

0.061 

0.310 

0,072 

0.139 

0,1100 

0.136a 

0.090 

0.157 

0.051 

0.093 

0.100 

0.814 


p.0U3 

-O.OQTf 

-O.pkfi 

-0.056 

-0.002 

-0.00 


-0.080 

-4.006 

-0.021 

0.0138 

-0,0835 

-0.013 

-0.015 

0.08$ 

0.016 

-0.036 

-0.091 

CW 

-0.112$ 

-0.1135 

-o.flol 

-o.ace 

-0.133 

-0.157 

-0.136 

-0.130 

-O.058 

-0.077 

-0.3299 

-0.3105 

-0.188 

-0.183 

-O.O67 

-0.069 

-0.318 

-0.389 

«Tp 

-0,06Efl 

-0.0361 

- - - 


-0.016 

0.001 

-0.016 

0.009 

-o.ool 

O.OOB 

-0.016a 

0.1150 


— 

- - - 

- - - 

- - - 

- - - 

°Ir 

Q.86S3 

0.2660 

- *• - 

» " *■ 

0.376 

0.376 

0.373 

0.373 

0.188 

0.1*4. 

0.1970 

0.1830 

- - - 

— 

- - - 

- - - 


— 


-0,3300 

-0.3370 

-0.690 

-0.690 

-0.619 

-0.619 

-0.638 

-0.630 

-0.318 

-0.362 

-0.5730 

-0.5730 

-0.673 

-0.690 

-0.590 

-0-595, 

-0.810 

-0,765 


■0.1660 

0.1670 

0.116 

0.116 

0.121 

0,121 

0.101 

0.109 

0.063 

0.069 

0,1011 

0.1866 

0.103 

0.110 

0.114 

0.110 

. 0.193 

0.171 

—5. 

-O.OT33 

■o.crra 

-0,107 

-0,105 

-o.no 

-0.13X5 

-O.C776 

-0,070 

-0,010 

-o.ow 

-O.06U 

-0.0726 

-0.064 

-0,102 

-0.02! 

-0.03 4 

-onto 

-0.210 
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Figure i -Two-view silhouettes of the nine airplanes investigated. 



Azimuth reference 



Figure 2. -System of axes and angles used in calculations. Arrows 
indicate positive values of forces, moments, and angles. 
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Calculated 

o / 6,300 "I 

□ 17, 150' l 

O 20,250' f 

A 24500' J 


Flight test 


Filled symbols indicate controls known fixed 
Half -filled symbols indicate controls assumed fixed 
Open symbols indicate one or more controls moving 



Figure 4. - Variation of period and time to damp to half 
amplitude with Mach number for airplane B. 
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Calculated 

O 10 , 000 ' 1 

D 25,000' > Flight test 

O 30,000 ') 

Half- filled symbols indicate controls 
assumed fixed 



Mach number 


Figure 5. - Variation of period and time to damp to 
half amplitude with Mach number for airplane C. 
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— — Calculated 

o 8 , 000 ' 1 

O 15,000' l Flight 

<5 15,000', H% added fin area on vertical tain test 

□ 30,000' J 

Filled symbols indicate controls known fixed 
Half -filled symbols Indicate controls assumed fixed 
Open symbols indicate one or more controls moving 



Mach number 




Figure 7.- Variation of period and time to damp to half 
amplitude with Mach number for airplane E . 
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Calculated 


• 10 , 000 '' ] 

■ 25,000' > Flight test 

♦ 30,000‘ J 

Filled symbols indicate controls known fixed 


<3 Jb 


!«• 




■S 

I 


I 



Mach number 


Figure 8. - Variation of period and time to damp to half amplitude 

with Mach number for airplane F. 
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• 10 , 000 '} 

■ 35,000’ \ 


Flight test 


Filled symbols indicate controls known fixed 



Q I I I I I I I I I I I I I I I 

.3 .4 5 £ .7 .8 .9 t.O 


Mach number 


Figure 9 Variation of period and time to damp to half amplitude 
with Mach number for airplane G. 
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Half -filled symbols indicate controls assumed fixed 


Open symbols indicate one or more controls moving 



Mach number 


Figure lO.-Variation of period and time to damp to half amplitude 
with Mach number for airplane H. 
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* 



o IO.OOO'\ 

A 20, OOO' > Flight test 

O 30,000' J 

Calculated for 10,000', neglecting C„ 

p 

Half-filled symbols indicate controls assumed 
fixed 



Figure H.~ Variation of period and time to damp to half 
amplitude with Mach number for airplane I. 
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Filled symbols indicate controls known fixed 
Half -filled symbols indicate controls assumed fixed 
Open symbols indicate one or more controls moving 



0 24602460246 


Calculated period, P, sec 


Figure 12.- Comparison between experimental and calculated period 
for the nine airplanes investigated . 
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Filled symbols indicate controls known fixed 
Half- filled symbols indicate controls assumed fixed 
Open symbols indicate one or more controls moving 



Calculated damping, t/j , per sec 


Figure 13- Comparison between experimental and calculated damping 
for the nine airplanes investigated . 
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